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ABSTRACT 

Aims. To investigate magnetic reconnection rates during the coalescence of two current loops in the solar chromosphere, by altering the neutral- 
hydrogen to proton density ratio, ioniziation/recombination coefficients, collision frequency and relative helicity of the loops. 
Methods. 2.5D numerical simulations of the chromosphere were conducted using a newly developed two-fluid (ion-neutral) numerical code. 
Developed from the Artificial Wind scheme, the numerical code includes the efi'ects of ion-neutral collisions, ionization/recombination, ther- 
mal/resistive difi'usivity and collisional/resistive heating. 

Results. It was found that the rates of magnetic reconnection strongly depend on the neutral-hydrogen to proton density ratio; increasing the 
density ratio by a thousand-fold decreased the rate of magnetic reconnection by twenty-fold. This result implies that magnetic reconnection 
proceeds significantly faster in the upper chromosphere, where the density of ions (protons) and neutral-hydrogen is comparable, than in the 
lower chromosphere where the density of neutral-hydrogen is over a thousand times the ion density. The result also implies that jets associated 
with fast magnetic reconnection, occur preferentially in the upper chromosphere / lower corona. The inclusion of ionization/recombination, 
an important physics efi'ect in the chromosphere, increases the total reconnected magnetic flux, but does not alter the rate of magnetic recon- 
nection. Reductions in the ion-neutral collision frequency, result in small increases to the rates of magnetic reconnection. The relative helicity 
of the two current loops was not observed to have any significant eifect on the rates of magnetic reconnection. Comparisons of two-fluid and 
MHD (MagnetoHydroDynamic) simulations, show significant difi"erences in the measured rates of magnetic reconnection, particularly for the 
higher neutral density cases which represent the lower chromosphere. This demonstrates that MHD is not an appropriate model for simulating 
magnetic reconnection in the solar chromosphere. 

Conclusions. The magnetic reconnection rates of coalescing current loops are strongly affected by the inclusion of neutral-hydrogen particles. 
It is therefore essential that ion-neutral collisions are included in future analytical/numerical models of chromospheric magnetic reconnection. 

Key words. Plasmas, Methods: numerical. Sun: magnetic fields. Sun: chromosphere. Sun: photosphere 



1. Introduction 

Magnetic reconnection is an important physical process in 
almost all cosmological and laboratory plasmas, because it 
serves as a mechanism for converting stored magnetic field 
energy into heat and non-thermal plasma energy. Magnetic 
reconnection in fully ionized plasmas has been extensivel 
investigated us ing the resistive-MHD jPriest & ForbesI uOO' 
iBisk amp 2005) and the Hall-MHD approximations, as well as 
in collisionless plasmas. The study of collisionless magnetic 
recon nection is well summarized by both the GEM (IBtrn et al. 
200 lb and Newton (IBirn et al.ll2005h reconnection challenges. 
Where in each challenge, a variety of different numerical codes 
were used to simulate the coalescence of magnetic islands; a 
Harris-type equilibrium. 
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Magnetic reconnection in weakly ionized plasmas is im- 
portant in the so lar ph otosphere/chromosphere (Sakai 1996|; 
Bula noy & Sakail Il998t iLitvinenkol Il999l; iFurusawa & Sakai 
120001: ISakaietal.' '2006 ), as wel l as in the interstellar 



medium ( Dorman & Kuls rud] 19951 ; IZweibel & Bran denburg! 



1997t iHeitsch & Zweibell I2003h . Solar observations of spec 



troscopic line shifts and non-thermal broadening, confirm that 
the photosphere, chromosphere and lower corona of the Sun 
are all in a highly dynamical state. These indicate that mass 
supply and wave propagation are both important mechanisms 
for transferring energy from the photosphere/chromosphere 
to the corona and ultimately the solar wind. Studies con- 
ducted by SUMER ( Solar U ltraviolet Measurements of Emitted 
Radiation; Wilhelm lTai] ([l995)), part of the SOHO mis- 
sion, clarified that observed explosive event s are bi-directional 



jets generated by magn etic recon nection (Innes et alj 11997 



Innes & Tota l 999; Rous sev & Ga lsgaard 2002). Observations 



have also shown evidence of magnetic reconnection in the 
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chromosphere a nd photosphere dBellot Rubio & Beckl 12005 ; 



Mclntoshl l2007h . where the neutral-hydrogen to proton den- 
sity ratio changes f rom one to 10^. Recent observations by 



Shibata et alj ( l2007h using the SOT (Solar Optical Telescope) 



instrument onboard the Hinode satellite, show ubiquitous 
anemone jets in the chromospheres of active regions; strong 
indirect evidence for the existence of chromospheric magnetic 
reconnection. 

Previously, in an attempt to better understand t he observed 
transi ent phenomena in the solar chromosphere, ISakai et al 



(120061) developed a simulation code to describe the dynamics 
of two-fluid (ion-neutral) plasmas. The two fluids, are coupled 
through proton/neutral-h ydrogen collisions, as well as through 
ionization/recombination. lSakai et al.l ( 120061) were successfully 
able to simulate the coalescence of of two counter-helical cur- 
rent loops (magnetic flux tubes) in the upper chromosphere. 
They found that the dynamics of two-fluid plasmas are quite 
different from those described by the single-fluid MHD approx- 
imation. During the simulated coalescence of the two current 
loops, they also observed a number of dynamical effects, in- 
cluding; proton heating, the formation of bi-directional proton 
jets and slow bi-directional plasma flows 

In this study, we will attempt to extend the ISakai et al 



(12006 *) investigation into the magnetic reconnection of coalesc- 
ing current loops. To achieve this, we have developed a new 
two-fluid numerical code, which we use to investigate how the 
magnetic reconnection rates of the coalescing current loops is 
altered by changing the; (a) neutral-hydrogen to proton density 
ratio, (b) ionization/recombination coefficients, (c) ion-neutral 
collision frequency, and (d) relative helicities of the two current 
loops. The results will be applicable to both chromospheric and 
photospheric magnetic reconnection studies. 

In Sect. 2 we present the two-fluid equations used to de- 
scribe the dynamics of ion and neutral particles. In Sect. 3 we 
describe the numerical scheme used in our newly developed 
two-fluid code. In Sect. 4 we describe the simulation model 
and initial conditions used in this study. In Sect. 5 we present 
our simulation results and finally in Sect. 6 we discuss the con- 
clusions that can be drawn from this study. 

2. Two-fluid Equations 

In this section we present the basic equations used to describe 
two-fluid (ion-neutral) plasmas. In what follows, the plasma 
density, pressure, velocity and magnetic field are given respec- 
tively by p, P, V and B, where the subscripts p and n refer to 
the ion (proton) and neutral fluids. We begin with the equations 
for the neutral fluid 

dp„ 



dt 



+ v-(p„y„) 



+ V ■ ipnVnVn) + VP„ 



d(pnVn) 

dt 

de 

-^ + ^-[{s„+P„)V„]-q„ = -53, 

ot 
where the neutral plasma energy is given by 

Pn ^ PnWn? 



r-1 



(1) 

(2) 

(3) 

(4) 



and the adiabatic constant is y = 5/3. The coefficients a,-, a^ 
and Q-f. refer respectively to the effects of ionization, recombi- 
nation and ion-neutral collisions. The neutral heat flux is given 
by qn - ^' (APn/Pn), whcrc A is the heat transfer constant. 
The source terms S i, S2 and S 3 seen above, collectively de- 
scribe the effects of ionization/recombination, ion-neutral drag 
and collisional heating 

S 1 = -Pp [arPp - aipnj , (5) 

^2 = a'cPpp„{y„-Vp^-pp(arPpVp-aip„V„y (6) 

S3 = acPpPn {Vn - F,,) Vp. (7) 
Next we present the equations for the ion fluid 



^+V-{ppVp) = 5:, 
-^^^ + V ■ {ppVpVp) -JxB + VPp = S2, 



(8) 



(9) 



dB 



— - V X (y„ X b) = /yV^B, (10) 



dt 



ds 



-^+'^-[{^p+pp)yp]-ip -s,, 



where the plasma energy is given by 



Pp PpWp? \Bi 
Sn = r + : + 



r-1 



2;Uo' 



(11) 



(12) 



the lorentz force by 7 = (V x B)/fio and 77 is the magnetic dif- 
fusivity. The ion heat flux, which includes thermal conduction 
and Joule heating, is given by qp = V^(APp/pp) + fiorjUl^. 

3. Numerical Scheme 

We use a newly developed two-fluid numerical code to sim- 
ulate the three-dimensional dynamics of the ion and neutral 
fluids present in the solar photosphere and chromosphere. In 
this numerical code, which we refer to as TwoYama, the fluids 
are coupled predominantly through ion-neutral collisions but 
also through ionization/recombination effects. TwoYama was 
developed using the pre viously proposed Artifi cial Wind (AW) 
numerical schem e (Sokolov etalJI 19991120021) . which was suc- 
cessfully used by ISakai et al.l (12006 ) to simulate the magnetic 
reconnection of coalescing chromospheric current loops. The 
AW scheme is based on the fundamental physical invariance, 
Galilean (or more generally Lorentz), of the governing plasma 
equations. It works by choosing the frame of reference such 
that the flow under consideration is always supersonic, thus 
trivializing the upwinding process and enabling highly simpli- 
fied forms of shock-capturing numerical schemes to be used. In 
practise this is achieved by adding an artificial velocity (hence 
Artificial Wind) to the flow under consideration. The result is a 
highly accurate shock-capturing numerical code, that can sim- 
ulate chromospheric problems substantially faster than other 
well known shock-capturing codes. 

The TwoYama code was designed specifically to make use 
of MPI (Message Passing Interface) parallelization, so that 
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large 3D simulations could be quickly simulated. The ion- 
ization/recombination source terms seen in the previous sec- 
tion, are solved using a time-splitting method at each half- 
time step, whilst the Joule heating, heat conduction, and mag- 
netic diffusion terms are solved explicitly inside the AW 
scheme. To ensure unconditional stability of the stiff ion- 
neutral drag t erms, we use a tim e-implicit di fferencing method, 
adapted from lStond d 1997b and lTothl(ll995h to include ioniza- 
tion/recombination. The update of the ion/neutral momentum 
after each timestep, Af, is given by 



R = 



Pp 



(^aiAtp„V„ - a,.AtppVp) + acAt[ppp„V„ -p„ppVj,) 



1 + Pp (OrAt + UiAt) + UcAt [pn + Pp) 



p";'v";'^PpVp + R 



«+lyn+l 



p • p — I^P ' P • ■"' A'n 'n —PnVn "■ (Ij) 

Following iTothI (119941) . we also use the implicitly solved mo- 
mentum change, R, to update the collisional heating source 
terms in the ion/neutral energy equations 



n+l 



= e"n + V„ ■ R, 



n+l 



= 4' -V„ R. 



(14) 



Continuous (zero-gradient) boundary conditions are used for 
the simulation boundaries. TwoYama solves the normalized 
forms of the equations seen in the previous section, according 
to the normalization given in the table below (nip is the proton 
mass). 

Table 1. Normalization Constants 



Symbol 


Quantity 


Constant 


Value 


t 


Time 


To 


= 1 s 


n 


Number density 


"0 


= 5 X 10'« m-^ 


T 


Temperature 


To 


= 8400 K 


P 


Density 


po = riQinp 


«8.4x 10-" kgm-3 


P 


Pressure 


Po = nokeTo 


« 5.8 X 10-3 Pa 


V 


Velocity 


Co = (Po/Po)"^ 


«8.3kms-' 


^, >', z 


Length 


A = Co To 


« 8.3 km 


B 


Magnetic field 


Bo = UlioPo)"^ 


« 1.2 X 10-* T 



4. Simulation IVIodel 



of radius a - 140A a; 1.5 Mm, are located parallel to the z-axis 
and are each assumed to be initially in independent equilibrium 
states (VP, - J i X B,). The initial magnetic field components 
of each loop (see Fig. 1) are given by 



B.i = qiB, {y - y„) la. By = -qiB-^ {x - x^) I a, 
B, = hiBme-'^''''"^', 



(15) 



where r; - [{x - x„)^ + (y - ycd^V^^, and the centers of the 
two current loops are (xci,yc\) = (350A, 512A) and {Xc2,yc2) - 
(674A, 5 12A). The initial ion pressure of each loop is given by 



Pm + Pm\^ 



2 2 



1 



,-2(r,laf 



while the initial ion density is similarly given by 



Ppi = Poo + Poo "7 - 



2 2 



- 1 



-2{r,laf 



(16) 



(17) 



The twist parameter, ^,, is fixed at q\ - q2 - 10. The he- 
licity parameter, hi - +1.0, determines the orientation of the 
magnetic field's z-component, B~. This parameter is important, 
as two types of magnetic reconnection can occur during the 
coalescence of two current loops; complete magnetic recon- 
nection (which we call the counter-helicity case) where the B^ 
components are anti-parallel, and partial magnetic reconnec- 
tion (which we call the co-helicity case) where the B^ compo- 
nents are parallel. The initial neutral density and pressure are 
set isotropically to a multiple of the background density, poo, 
according the specific neutral-hydrogen to proton density ratio 
used in each simulation. The initial velocity of both the ion and 
neutral fluids is zero. 

The heat transfer constant, seen in the heat flux terms 
of Sect. 2 is held fixed in all simulations at /I - 1.5 x 
IQ-^ kg m-' s-'. The thermal conductivity coefficient can be 
obtained from the heat transfer constant from k - ksA/mp w 
0.012 W m-' K-' . The ion-neutral collision frequency is given 



by v^ 



pn 



^ p^ pn^tp: 



where the collision cross-section, cr, 



5.5 X 10 m , and the ion thermal velocity, v,p 



pn 



10.8 km s . Therefore 



300 Hz, giving a normalized 



ion-neutral collision coefficient Uc - OcPo^q - 300, where 
ffc = Vpnipp. If we assume that radiative recombination is dom - 
inant, then using the formula by Kap lan & Pikelnen (119701) : 



In this study, we assume a plasma temperature oiTp = T„ = a^. = 1.6 x 10^(10'*/r)'^-^'' m^ kg-' s-', which gives a nor- 
malized recombination coefficient of a^ - o'rPoTo ~ 0.015. 
Initially, we assume that recombination balances with ioniza- 
tion, Or = 5; = aipoTo and that pp - p„. 

Reconnection of magnetic fields occurs due to a non-zero 
magnetic diffusivity (or a finite conductivity). If we take the 
photospheric el ectron-hydrogen co llision frequency as v^,, w 
lO** - 10'" Hz dVranies et al.ll2008b . then we can estimate the 
magnetic diffusivity due to colUsions between electrons and hy- 
drogen atoms as 



8400 K and a plasma number density of Hp = n„ = 5x10 m"^. 
Using these parameters, the background sound speed Cs = 
iyPpIpp^^^ ~ 10.8 km s-'. For simplicity, we choose a nor- 
malization time of To = 1 s, such that the grid spacing A ^ 
10.8 km. Therefore the system resolution of A^j^ - Ny = 1024 
corresponds to a system size of L^ - ANx - Ly - ANy ^ 
11.0 Mm. The initial background magnetic field, pressure and 
density are given respectively by Boo = 1.0 x 10-"* T, Pqo = 
HphsTp K 5.8 X lO-"* Pa and poo = ripWip k 8.4 x 10-" kg m-^. 
Therefore, the corresponding background Alfven and fast- 
magnetosonic speeds are respectively c„ - Bo/ifioppy^^ ^ 
9.8 km s-i and c,„ = (c/ + Ca^'^ ^ 14.5 km s-'. 

We simulate the 2.5D coalescence of two current loops in 
the upper chromosphere. The two current loops (/ - 1,2), each 






5.6x10 



5-6 



Ks-'). 



(18) 



If we take the characteristic length scale as A ^ 10.8 km, 
and the characteristic velocity as the sound speed Cj w 
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10.8 km s ', then the magnetic diffusivity corresponds to a 
magnetic Reynolds number R,„ - Acjr] «; 20 - 200. In this 
study 7/ = 1.0 X 10^ m^ s"', which corresponds to /?,„ ^10. 

To gauge the rate of magnetic reconnection during the coa- 
lescence of the two current loops, we measure the total recon- 
nected magnetic flux, i/r, at the point x - y - 512A, using the 
following expression 



■JE.. 



(19) 



1//= - \ E,dt, 

where £. is the induced electric field. 

5. Simulation Results 



In this section we present the results of our numerical simula- 
tions. Unless stated, all of the results are given in normalized 
units, according the normalization specified in Table 1 . 

We begin with Fig. [1] where we show the initial magnetic 
field configuration for two co-helical current loops. The initial 
ion pressure in-balance between the two loops (Fig.|2]), causes 
them to spontaneously collapse towards each other. During the 
initial coalescence phase (t < 1000), strong inflows (Fig. [3]) 
develop along the x-axis, with subsequent outflows along the 
y-axis. This inflow of plasma compresses oppositely directed 
magnetic field lines, seen in the central region of Fig. |4] to 
form a strong reverse current sheet (see Fig. |5]) parallel to the 
y-axis. The current sheet is then dissipated by the plasma's fi- 
nite conductivity, leading to reconnection of the (B^ - By) mag- 
netic field lines. The significant heating (Joule and collisional) 
generated during the reconnection process, is clearly visible in 
Fig. |6l where we see a 30% rise in ion temperature along the 
current sheet. Conversely, in Fig.|2]we see a 5-10% reduction in 
neutral temperature along the current sheet, due to collisional 
cooling. This implies that 20-25% of the increase in ion tem- 
perature is due to Joule heating, while a further 5-10% increase 
is due to collisional heating with neutral particles. In Fig. [8] 
we see that during the initial coalescence phase, reconnection 
of magnetic fields lines reduces the magnetic energy, or mag- 
netic pressure, relative to the ion thermal pressure. This is turn 
accelerates the ion and neutral fluids (due to ion-neutral colli- 
sions) into the central region; an effect clearly seen in Fig. [8] 
by the rapid increase in ion and neutral kinetic energy. In the 
later coalescence stage {t > 1000), collisional heating is seen 
to slowly convert the ion and neutral kinetic energy into inter- 
nal (thermal) plasma energy. Fig.[8]demonstrates that the main 
effect of the magnetic reconnection process is the conversion 
of magnetic energy into kinetic and thermal energy. 

5.1. Effects of ion-neutral density ratio 

Firstly, we investigate how the rate of magnetic reconnection 
is altered by increasing the neutral-hydrogen to proton den- 
sity ratio, Pnlpp. We present numerical results for the coales- 
cence of two co-helical current loops, using density ratios of 
PnIPp - 1, 10, 100, 1000. The normalized collision frequency 
is held fixed at Oc - 300, while ionization and recombination 
effects are ignored (a,- - Ur - 0.0). 



It is well known that the neutral-hydrogen to proton density 
ratio increases with depth in the chromosphere/photosphere. 
Changing this density ratio therefore effectively corresponds 
to a change in depth at which the simulated current loop co- 
alescence, and subsequent magnetic reconnection takes place. 
The density ratio p„ jpp = 1 is representative of the upper chro- 
mosphere, while PnIpp - 1000 is representative of the lower 
chromosphere. In Fig. |9l we use a normalized form of Eq. 19 
to plot the total reconnected magnetic flux as a function of time, 
for various neutral-hydrogen to proton density ratios. The slope 
of each line is a direct measure of the magnetic reconnection 
rate (as measured ?Lt x - y - 512A). From this figure, it is im- 
mediately clear that increasing the neutral-hydrogen to proton 
density ratio strongly decreases the magnetic reconnection rate. 
Comparing the lines representing p„/pp = 1 and p„/pp = 10, 
we see that a ten-fold increase in density ratio reduces the mag- 
netic reconnection rate by a factor of two. Whilst, comparing 
the lines representing p„ /pp = 1 and p„/pp = 1000, we see that 
a thousand-fold increase in density ratio reduces the magnetic 
reconnection rate by a factor of twenty. The simulation results 
seen in Fig. |9] therefore suggest that, in the case of coalescing 
current loops, magnetic reconnection proceeds faster in the up- 
per chromosphere than in the lower chromosphere. 

Fig. |9] shows a direct correlation between magnetic recon- 
nection rate and density ratio, however the correlation between 
the peak reconnected magnetic flux, i/r, and the density ratio is 
less clear. The peak value for the lines p„/pp = 1 and Pn/Pp - 
10 are approximately equivalent at i/' « 21, whereas the lines 
representing the higher density ratios of p„/pp = 100 - 1000, 
appear to indicate a reduction in the peak reconnected magnetic 
flux with increasing density ratio. 

Finally, comparisons of magnetic and velocity vector plots 
for various density ratios, in addition to Fig.|9l appear to show 
that the density of neutral particles does not alter the funda- 
mental dynamics of coalescence process itself, but simply the 
rate at which it occurs. 



5.2. Effects of ionization/recombination 

The recombination of protons and electrons to form neutral 
hydrogen, is an important physical process occurring in the 
partially ionized upper chromosphere. We therefore investigate 
how the magnetic reconnection rate is altered by including the 
effects of ionization and recombination. We show results for 
two co-helical current loops using a density ratio of p„/pp = 1. 
The normalized collision frequency is fixed at a^. = 300, while 
the normalized ionization and recombination coeflicients are 
set to a, = 5, = 1.5,0.15,0.015,0.0. 

In Fig. [To] for clarity we do not show the two cases of 
&! - Or - 0.15,1.5 since these gave visually indistinguish- 
able results to the a, - a^ - 0.015 case. From the figure we 
see that the inclusion of ionization/recombination effects has 
no effect on the rate of magnetic reconnection. We do however 
see an increase in the total reconnected magne tic flux by ^ 5%. 



This phenomena was previously explained by ISakail (Il996h as 
resulting from a decrease in ion pressure at the current sheet, 
formed during the coalescence process (see Fig. |5]i. This de- 
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crease in ion pressure results from protons and electrons in the 
current sheet recombining to form neutral hydrogen. Since it is 
the magnetic field which drives the reconnection process, any 
reduction of the ion pressure, which alone balances the mag- 
netic pressure, will result in enhanced magnetic reconnection. 

5.3. Effects of collision frequency 

Next we investigate how the magnetic reconnection rate is al- 
tered by changing the ion-neutral collision frequency. We show 
simulation results for two co-helical current loops using nor- 
malized collision frequencies of 5-^ = 300,3.0,0.03,0.0. The 
density ratio is fixed at Pn/Pp = 1, while ionization and recom- 
bination effects are ignored (a,- = a,- = 0.0). 

Firstly, in Fig.[TTl we note that the lines representing Uc = 
300 and o-e = 3 produce nearly identical results (only ac = 300 
is visible). This is because both cases represent strongly cou- 
pled ion-neutral plasmas; the ion-neutral collisions cause the 
ion and neutral particles to behave as a single heavy fluid with 
a combined ion-neutral density. For the a^- = 0.0 case on the 
other hand, the ion and neutral particle dynamics are com- 
pletely decoupled. The ions therefore move as a single light 
fluid (equivalent to a MHD simulation). The collapse and even- 
tual coalescence of the two current loops is driven by a continu- 
ous reduction in magnetic pressure (due to magnetic reconnec- 
tion), balanced only by the ion thermal pressure. The neutral 
particles by definition are not influenced by the magnetic field, 
and therefore play no part in balancing the magnetic pressure. 
However, the eff'ective inertia of the initially stationary fluids, 
does depend on both the neutral-hydrogen to proton density 
ratio and the ion-neutral collision frequency. Given that the ini- 
tial conditions for the ion particles and magnetic field are the 
same in each simulation, the initial acceleration of the fluids 
can only therefore depend on their eff'ective inertia. In strongly 
coupled ion-neutral plasmas, where the eff'ective inertia is rel- 
atively high, we would expect the coalescence and subsequent 
magnetic reconnection to proceed, at least initially, relatively 
slowly. Whereas in weakly coupled ion-neutral plasmas, where 
the effective inertia is relatively low, we would expect the re- 
verse to be true. This appears to be confirmed by Fig. (TT] where 
we see a slight increase in the rate of magnetic reconnection, 
with decreasing ion-neutral collision frequency. 

In Fig. [n] the rate of magnetic reconnection for the MHD 
case (ffc = 0.0), does not diff'er significantly from the two-fluid 
cases (q-c = 300, 3.0, 0.03). We note however, that this is only 
true for the low neutral density case, where p„/pp - 1. If we 
compare the MHD c ase of Fig . [TT] to the p„ /pp - 1 000 c ase of 
Fig. |9] we see a twenty-fold difference in the magnetic recon- 
nection rates. The MHD approximation is therefore not an ap- 
propriate model to describe magnetic reconnection in the lower 
chromosphere or photosphere. 

Finally in Fig. [TT] we note that the peak reconnected mag- 
netic flux appears to decrease with decreasing ion-neutral colli- 
sion frequency. We speculate that the increased effective inertia 
of ion-neutral fluids may extend the time for loop coalescence, 
thus leading to an increased reconnected magnetic flux. This 
effect is of course balanced by the reduced initial acceleration. 



caused by increased neutral density, discussed above. The re- 
sults of Fig. [9] would therefore suggest that neutral-hydrogen to 
proton density ratios of p„/pp = 1-10, would lead to optimum 
amounts of reconnected magnetic flux. 



5.4. Effects of loop helicity 

In this final subsection, we investigate how the magnetic re- 
connection rate is altered by changing the relative helicities 
of the two current loops. We show simulation results for two 
counter-helical (opposite B,) current loops, using density ra- 
tios of p„/pp - 1, 10, 100, 1000. The normalized collision fre- 
quency is held fixed at a,. = 300, while ionization and recom- 
bination effects are ignored (or,- - Ur = 0.0). 

In Fig. [12] we see that, as previously seen in the co-helical 
case, the rate of magnetic reconnection strongly depends on 
the density ratio. Comparisons of Fig.[T2]to the co-helical case 
seen in Fig. [9] show only minor variations in the both the rates 
of magnetic reconnection and the total reconnected magnetic 
flux. We expect that this result holds true only for the 2.5D case 
studied here, since previous 3D studies have shown that the 
relative helicities of the loops does indeed affect the magnetic 
reconnection process. A subsequent 3D study will investigate 
the relative helicity further 

6. Discussions and Conclusions 

In this study, we investigated the magnetic reconnection of 
two coalescing chromospheric current loops. We conducted 
2.5D simulations of the coalescence using a newly developed 
two-fluid (ion-neutral) code, which includes the effects of ion- 
neutral collisions, ionization/recombination, thermal/resistive 
diffusivity and coUisional/resistive heating. 

In the results of the previous section, we found that the fun- 
damental physical dynamics of the coalescence process was not 
altered by changing (a) the neutral-hydrogen to proton density 
ratio, (b) the ionization/recombination coefficients, (c) the ion- 
neutral collision frequency or (d) the relative helicities of the 
two loops. We did find however, that the rate of magnetic re- 
connection was susceptible to each parameter to a varying de- 
gree: 

(a) Increasing the neutral-hydrogen to proton density ratio 
strongly decreased the rate of magnetic reconnection. 

(b) The inclusion of ionization and recombination effects did 
not alter the rate of magnetic reconnection. 

(c) Decreasing the ion-neutral collision frequency slightly in- 
creased the rate of magnetic reconnection. 

(d) The relative helicities of the two current loops had no dis- 
cernable effects on the rate of magnetic reconnection. 

Altering the neutral-hydrogen to proton density ratio, had a sig- 
nificantly stronger effect on the rates of magnetic reconnection, 
than any of the other parameters listed above. In fact, we found 
that the rate of magnetic reconnection for p„ jpp - 1 (upper 
chromosphere) was twenty times faster than for p„/pp = 1000 
(lower chromosphere). This important result, which may be 
applicable to other forms of chromospheric magnetic recon- 
nection, implies that observed jets associated with fast mag- 
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netic reconnection, occur preferentially in the upper chromo- 
sphere / lower corona. Comparisons of two-fluid and MHD 
(MagnetoHydroDynamic) simulation results, show significant 
differences in the measured rates of magnetic reconnection, 
particularly for the density ratios representative of the lower 
chromosphere. This demonstrates that MHD is not an appropri- 
ate model to describe magnetic reconnection in the solar chro- 
mosphere/photosphere. 

The main result of this study shows that the rates of mag- 
netic reconnection for coalescing current loops, are strongly 
afifected by the inclusion of ion-neutral collisions. We therefore 
conclude that it is essential that the effects of neutral particles 
are included in future analytical/numerical models of chromo- 
spheric magnetic reconnection. 
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Fig. 1. The initial magnetic field: a contour plot of B^ overlayed 
by a vector (B^ - By) plot- 
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Fig. 2. Time development of the ion pressure for two co-helical 
current loops, where p„/pp = 1. 
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Fig. 3. Time development of the ion velocity (Vx - Vy) for two 
co-helical current loops, where p„/py, - 1. 
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Fig. 4. Time development of the magnetic field (B^ - By) for Fig. 7. Time development of the neutral temperature for two 



two co-helical current loops, where pn/pp = 1. 



co-helicalcurrent loops, where p„/pp = 1. 
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Fig. 5. Time development of the current density, J^, for two co- Fig. 8. Time development of plasma energy for two co-helical 



helical current loops, where p„/pp = 1. 
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current loops (pn/Pp = 1)- Kinetic and magnetic energy plots: 
(solid) x-component, (dotted) y-component and (dashed) z- 
component. Internal energy plot: (solid) ions and (dashed) neu- 
tral. Each energy component is normalized to the total plasma 
energy. 
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Fig. 6. Time development of the ion temperature for two co- 
helical current loops, where p„/pp = 1. 
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Fig. 9. Effect of neutral-hydrogen to proton density ratio on Fig. 11. Effect of ion-neutral collision frequency on the recon- 

the reconnected magnetic flux for two co-helical current loops, nected magnetic flux for two co-helical current loops (pn/Pp - 

(soUd) p„/pp - 1, (dotted) p„/pp - 10, (dashed) pn/Pp - 100 1). (soHd) Oc - 300, (dotted) Uc - 3.0, (dashed) Oc - 0.03 and 

and (dot-dashed) p„/pp = 1000. (dot-dashed) a^ = 0.0 (MHD). 




Fig. 10. Effect of ioniziation/recombination on the the recon- Fig. 12. Effect of neutral-hydrogen to proton density ratio on 
nected magnetic flux for two co-helical current loops (p„/pp = the reconnected magnetic flux for two counter-helical current 
1). (solid) Oi - ar - 0.0 and (dotted) a, - Ur - 0.015. loops, (solid) p„/pp - 1, (dotted) p„ /p^ = 10, (dashed) p„/pp = 

100 and (dot-dashed) p„/pp = 1000. 



